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Recent progress on the event-by-event analysis of intermittent data by R. A. Janik 
and myself is reported. The intermittency analysis of single event data (particle 
moments) in multiparticle production is improved, taking into account corrections 
due to the reconstruction of history of a particle cascade. This approach is tested 
within the framework of the a-model. 



1 Introduction 



The first data on possibilc intermittent behaviour in multiparticle production 
EJ came from the analysis .of the single event of high multiplicity recorded by 
the JACEE collaboration □. It was soon realized, however, that the idea may 
be applied to events of any multiplicity provided that averaging of the distri- 
|| butions is performed H. This led to many successful experimental studies of 

intermittency □, and allqwed to express the effect in terms of the multiparti- 
cle correlation functions It should be realized, however, that the averaging 
procedure, apart from clear advantages, brings also a danger of overlooking 
some interesting effects if they are present only in a part of events produced 
in high-energy collisions. For example, the unique properties due to the pres- 
ence of quark-gluon plasma in multiparticle production would manifest only in 
some events, see e.g. 0. Taking into account the sample of events and averag- 
ing, oyer them destroys such an information. Therefore, as already discussed 
inp, □, there is a need for event-by-event analysis of multiparticle production 
processes. In this way the fluctuations of the measured physical quantities (e. 
g. factorial moments) from event to event can be observed and estimated, and 
any anomalous behaviour of them has a chance to manifest very clearly. Such 
studies should necessarily be restricted to high- multiplicity events because only 
there one may expect the statistical fluctuations to be under control. 

SucIl aiyapprpach to the multiparticle data analysis has been already pro- 
posed in0,l3,lj. InLj a new quantity: erraticity has been introduced to investigate 
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the event-by-event fluctuations of factorial moments, and to search for their 
properties. Erraticity denotes the normalized moment of event-by-event dis- 
tribution of a horizontally averaged factorial moment. It probes both types 
of fluctuations: horizontal ones connected with the spatial bin pattern and 
vertical ones i.e. event-by-event ones. 

Ini the event-by-event fluctations of particle moments have been investi- 
gated directly for the one-dimensional a— model of random cascading. Monte- 
Carlo simulations of the model allowed one to obtain the histograms of event- 
by-event distributions of horizontally averaged particle momenta and estimate 
the relation between the intermittency parameters obtained from such a his- 
togram, and the intermittency parameters derived after usual procedure of 
averaging particle moments over all events. The results were promising: the 
average value of the intermittency exponent reproduced well the value obtained 
by averaging particle moments over events, however with the tendency to un- 
derestimate the theoretical value. Furthermore, the dispersion of the moment 
distribution was inversely proportional to the length of a generated cascade, 
and even for short cascades substantially smaller than the average value. The 
latter property was of a special importance : it allowed one to distinguish be- 
tween groups of events emerging from cascades with different characteristics. 

In this paper we would like to improve the analysis of single event data 
presented inlj. Taking into account corrections due to the method of recovering 
the history of the multiparticle cascadeu,l3, we expect to reduce the discrepancy 
between the theoretical value of intermittency exponent and its value estimated 
from the event-by-event histogram □. Our discussion will proceed as follows. 
In section 2 we recall the definition of the intermittency exponents and the 
technique used to calculate them El, □. In section 3 the definition of the a model 
will be briefly presented, and applied in section 4 to calculate corrections for 
extracting intermittency exponents from single event data. Section 5 is devoted 
to the comparison of theoretical results with numerical simulations. Finally in 
section 6 we present our conclusions. 

2 Intermittency exponents 

Consider a multiparticle production cascade distributed into M bins. At the 
n th stage of the cascade we measure the distribution of a particle density into 
M bins. Assume for simplicity that M = 2". We thus have 2™ numbers 
(quantities) denoting the content of each bin : 



(1) 



To perform the event-by-event analysis one is interested in the behaviour of 
particle moments with the stage of the cascade : 

i=0 

The scaling_hehaviour of these moments is parametrized by intermittency ex- 
ponents g El : 

4«)~2 n ^«. (3) 

The task is to estimate the value of an intermittency exponent. There are two 
different ways of doing it. The first one is to calculate the average moment z\ 
for the whole ensemble of individual events, and from this to reconstruct the 
intermittency exponent. The second one is to calculate the exponent (\> q for 
each event separately, and then to recover the average 4> q . The latter approach 
has the advantage of being able to distinguish between two independent cas- 
cading processes each with different <fi q . This could be done by looking at the 
distribution of individual <fi q 's. In the former method both of these possibly 
independent processes would be artificially forced to be described by a single 
'effective' 4> q . 

In the following we would like to address the question of reliably recon- 
structing the correct value of 4> q from single event data. Numerical simulations 
inEI showed that there is an inherent discrepancy between the theoretical value 
and the distributions of event-by-event 4> q (see Tables 1,2). The aim of this 
letter is to analyze this result and introduce a correction which improves the 
estimation. 

A convenient way of calculating <j) q is to make a linear fit to the points 
(n,log4™^) ( au logarithms are taken to be calculated in base 2, i.e. logx = 
lnx/ln2) : 

log 4") =n-4> q +b. (4) 

This procedure has the advantage of cancelling out the major part of the correc- 
tion coming from the fact that we are effectively reconstructing the exponents 
from (logZg™^) while the true value is defined in terms of log(zq™^). 

However there is still one caveat to (|J) . Since we cannot in general separate 
out the various stages of the cascade, one reconstructs the previous stages from 
the last one by summing the x\ s in adjacent-bins using the technique de- 
scribed inty ( and applied there to JACEE event □ ). Namely one approximates 
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the true value of xf 1 fe ' by : 



~ 2 k ^ X 2*y.i+j- ( 5 ) 
3=0 

Therefore in (m one really uses the reconstructed moments : 



(n—k) (n—k) 

A — > v\ = 2 , 

3=0 



jE(»r')' w 



constructed 

3=0 

We will now use the a model of random cascading to calculate explicitly 
the difference between the true and reconstructed moments and the resulting 
shift of the <j) q distribution from the theoretical value. 

3 The a model of random cascading 

In the a model of random cascading§ the root of the cascade — x^ is taken 
to be a with probability p a and b otherwise (with probability pb = 1 — p a ). 
One generates the next stages of the cascade recursively. The two bins x^ +1 ^ 
and x^Y^ are obtained from by : 

I with probability p a , (7) 

a4i — ► b ■ x[ n ^ with probability pb, (8) 

and same for x^t^ . The parameters a and b are taken to satisfy : 

ap a + bpb = l. (9) 

Particle moments fulfill the relation : 

z (n) = 2 (n+i)-0 < , ! ( 1Q ) 

where intcrmittency exponents <p q are equal to : 

<^=log(a^ a + &<%). (11) 

4 Reconstructed moments 

The reconstructed moments in the a model are related to the true ones by : 



9 

' 'q, reconstructed ^r. 



d = ^ E ( )^z q n - k) - Pq (k). (12) 



i=0 \ j=0 
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where the average (...) is taken over the random choices made only above the 
(n — fc)-th stage of the cascade. The factor p q (k) can be calculated exactly 
(see below) and we propose to use it to compensate for the errors introduced 
by the reconstruction procedure. In particular the reconstructed moments 
entering will be shifted by : 

1°S Z q, reconstructed * Z q, reconstructed ~ l°g(Pg(^))- (13) 

We will now determine the explicit form of the correction p q (k). By the 
definition of the a model, the correction p q (k) can be calculated just by eval- 
uating : 

p«( k ) = (Uk X>fj ) ( 14 ) 

in the a model modified by taking the starting bin = 1. 

First it is easy to see that for q = 1 there is no correction pi(k) = 1. This 
is due to (^). Also all corrections vanish for k = : 

p,(0) = 1. (15) 

The appearance of a correction for q > 1 comes from the fact that the 'number' 
of particles in this model has a nonzero dispersion. 

Consider first the case of q = 2. We will now split the bins (xi's) appearing 
in Q into a left half (i < 2 k - 1 ) and a right half (i > 2 k ~ 1 ): 




+2(iE' t )(^))- (16) 



Using the fact that the left and right bins are independent one gets the recur- 
rence relation: 

Pa(*) = I ( Pa a 2 +p b b 2 )p 2 (k - 1) + i (17) 



This can be solved together with the initial data ([0]), to yield a closed form 
solution: 

/ x f d 2 \ k 1 -do 1 
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In general one can obtain the recurrence relation for general q in exactly 
the same way: 



*<*>-5i£(!) 



did q -i -pi(k - l)p q -i(k - 1) 



(19) 



where 



di = p a a l +Pbb l . 



(20) 



5 Discussion 

We have performed numerical simulations of the a-model in order to test the 
improved single data analysis in practice. In Fig. 1 and Fig. 2 the histograms 
of the corrected (with the shift (|l^) taken into account) and standard (without 
the correction (p^)) values of intermittency exponents tp2)<P3 are plotted for 
90000 generated cascades of 5 and 10 steps. The peaks with the correction 
included are significantly closer to the theoretical value. The dispersion of the 
distribution estimated directly from the observed peak, for the "corrected" his- 
togram is smaller than the dispersion of the "standard" one. It decreases with 
the number of cascade steps. The numerical values of "corrected" and "stan- 
dard" dispersion as a function of the cascade length are presented in Tables 1, 
2 for 2 different sets of cascade parameters. The corrected dispersion is rela- 
tively small, and it allows to distinguish between the cascades with different 
parameters (Fig.l, 2). 

The influence of the correction (^|) on the value of the intermittency ex- 
ponents obtained from averaging over the ensemble of events ('center of mass' 
of the histogram) was also investigated. The results are presented in Tables 
3,4 for 2 different sets of cascade parameters. The estimation of intermittency 
exponents for the corrected case is much better than for the standard one. 

In the preceding, the formula for the correction (see e.g. (|l8|)) depends on 
the values of the parameters a, b of the a-model. In practice, however, one 
would like to implement some sort of model independent correction. A possible 
way of doing this is to use the fact that the corrections logp2(i) and \ogps(i) 
seem to change most dramatically in the first few steps of the reconstruction 
procedure (near the 'end' of the cascade). After that they seem to stabilize at 
some constant value. This would suggest using just the reconstructed moments 
near the beginning of the cascade in the fit (^). In practice, however, this might 
perhaps suffer from low statistics and large fluctuations. 

An alternative procedure would be to first determine the parameters a and 
b using the standard (uncorrected) method, and then substitute those param- 
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eters into ( |19| ) and use the improved analysis to obtain a better approximation 
of the exponents. One could repeat this until the result no longer changed. 

6 Conclusions 

Our conclusions can be summarized as follows : 

(a) the value of intermittency exponent estimated from the maximum of 
"corrected" histogram moves closer to the theoretical value, 

(b) the dispersion of the distribution estimated directly from the observed 
peak for the "corrected" histogram is smaller than the dispersion of the "stan- 
dard" one, 

(c) the corrected value of intermittency exponent obtained after averaging 
over the sample of events estimates the theoretical value better than in the 
standard case, 

(d) a possible procedure of improving the analysis without the knowledge 
of a-model parameters is proposed. 

Acknowledgments 

We would like to thank A. Bialas for discussions and suggestions. This work 
was supported by Polish Government grant Project (KBN) grant 2P03B00814 
and 2P03B04214. RAJ and BZ were supported by the Foundation for Polish 
Science (FNP). 

References 

1. A. Bialas, R. Peschanski, Nucl. Phys. B273, 703 (1986) 

2. JACEE coll. T. H. Burnett et al. , Phys. Rev. Lett. 50 (1983) 2062 

3. B. Buschbeck, P. Lipa, Phys. Lett. B223 (1989) 465 and private com- 
munication 

4. for recent reviews see 

E. A. de Wolf, I. M. Dremin, W. Kittel, HEN-362 (1993) update July 
1995, 

P. Bozek, M. Ploszajczak, R. Botet, Phys. Rep. 252 (1995) 101-176 

5. P. Carruthers, I. Sarcevic, Phys. Lett. B189 (1987) 442 

6. A. Bialas, R. Hwa, Phys. Lett. B253 (1991) 436 

7. A. Bialas, A. Szczerba, K. Zalewski, Z. Phys. C46 (1990) 163 

8. R. Hwa, Phys. Lett. B253 (1991) 436 

9. A. Bialas, B. Ziaja, Phys. Lett. B378 (1996) 319 



7 



Figure captions 



Fig. 1, 2 Histograms of the intermittency exponents fa (left column) and fa 
(right column) simulated for the set of parameters a = 0.8, b = 1.1 (Fig. 1) 
and a = 0.5, b — 1.5 (Fig[B. 2) in 90000 events for 5 and 10 cascade steps. The 
wider curves correspond to 5 stages of the cascade. 'Solid' curves represent the 
histograms with the correction (Il9|) taken into account. 
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Table 1. Standard and corrected intermittency exponents (determined from 
the position of the maximum of the histograms) and their dispersions (errors) 
for a = 0.8, b = 1.1 and n = 5, . . . , 10 cascade steps. Theoretical values for 
intermittency exponents are <p2,theor — 2.85 x 10~ 2 and ipz^heor = 8.13 x 10~ 2 . 



ipi = lQ-' A x 


5 


6 


7 


8 


9 


10 


(f2 


1.90 ±0.89 


2.00 ±0.82 


2.26 ±0.75 


2.49 ±0.66 


2.52 ±0.59 


2.46 ±0.52 


(P2.corr 


2.66 ±0.75 


2.79 ±0.66 


2.66 ±0.56 


2.85 ±0.46 


2.85 ±0.43 


2.85 ±0.36 


(f3 


5.66 ±2.46 


5.66 ±2.46 


6.48 ±2.05 


6.64 ± 1.72 


6.81 ±1.72 


6.72 ± 1.49 


(PZ.corr 


7.79 ±2.13 


7.63 ± 1.81 


7.62 ± 1.64 


8.00 ± 1.40 


8.12 ± 1.15 


8.12 ±0.98 



Table 2. Intermittency exponents (determined from the position of the maxi- 
mum of the histograms) and their dispersions (errors) for a = 0.5, b = 1.5 and 
n = 5, . . . , 10 cascade steps. Theoretical values for intermittency exponents 
are tp2,theor = 3.22 x 10 _1 and ip3,theor = 8.07 x 10 _1 . 



(Pi = 10 _i x 


5 


6 


7 


8 


9 


10 


(f2 


2.00 ± 1.00 


2.09 ±0.92 


2.43 ± 0.74 


2.. 61 ±0.70 


2.44 ±0.74 


2.26 ±0.70 


(P2.corr 


3.13 ±0.79 


3.13 ±0.74 


2.96 ±0.61 


2.87 ±0.57 


3.05 ±0.52 


3.00 ±0.48 


(fi3 


4.52 ±2.23 


5.31 ±2.23 


5.83 ±2.02 


5.90 ± 1.83 


6.03 ± 1.83 


5.70 ± 1.70 




7.53 ± 1.90 


7.93 ± 1.83 


7.66 ± 1.57 


7.53 ± 1.31 


7.66 ± 1.31 


7.66 ± 1.18 
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Table 3. Standard and corrected intermittency exponents and their dispersions 
(errors) for a = 0.8, b = 1.1 and n = 5, . . . , 10 cascade steps obtained after 
averaging over the sample of 90000 events. Theoretical values for intermittency 
exponents are (p2,theor — 2.85 x 10~ 2 and ip^jheor = 8.13 x 10~ 2 . 



ipi = lQ-' A x 


5 


6 


7 


8 


9 


10 




2.16 ±4.34 


2.45 ± 1.51 


2.57 ±0.74 


2.63 ±0.66 


2.68 ±0.60 


2.72 ±0.54 


^P2.corr 


2.90 ±0.70 


2.89 ±0.58 


2.88 ±0.49 


2.88 ±0.42 


2.87 ±0.36 


2.87 ±0.31 


^3 


6.54 ±4.80 


7.03 ±2.62 


7.33 ±2.11 


7.53 ± 1.90 


7.68 ± 1.73 


7.77 ± 1.58 


(fS.corr 


8.32 ±2.00 


8.29 ± 1.67 


8.25 ± 1.40 


8.23 ± 1.18 


8.22 ± 1.02 


8.20 ±0.89 



Table 4. Intermittency exponents and their dispersions (errors) for a = 0.5, b = 
1.5 and n = 5, . . . , 10 cascade steps obtained after averaging over the sample 
of 90000 events. Theoretical values for intermittency exponents are if2,theor = 
3.22 x 10" 1 and </? 3 , t/leor = 8.07 x 10" 1 . 



(Pi = 10 _i x 


5 


6 


7 


8 


9 


10 




2.33 ± 1.2 


2.50 ±0.95 


2.62 ±0.82 


2.69 ±0.76 


2.77 ±0.74 


2.81 ±0.72 


i P2.corr 


3.20 ±0.70 


3.17 ±0.63 


3.15 ±0.57 


3.15 ±0.52 


3.14 ±0.47 


3.14 ±0.43 


f3 


5.78 ±2.36 


6.13 ±2.06 


6.38 ± 1.90 


6.55 ± 1.81 


6.71 ±1.78 


6.81 ± 1.75 


tP3,corr 


8.16 ± 1.71 


8.05 ± 1.51 


7.96 ± 1.36 


7.90 ± 1.24 


7.86 ± 1.13 


7.84 ± 1.06 
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Figure 1: a-b Histograms of the intermittency exponents a) 02 and b) 03 
simulated for the set of parameters a = 0.8, b = 1.1 in 90000 events for 5 and 
10 cascade steps. The wider curves correspond to 5 stages of the cascade. 
'Solid' curves represent the histograms with the correction p q (k) taken into 
account. 
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Figure 2: a-b Histograms of the intermittency exponents a) 02 and b) 03 
simulated for the set of parameters a = 0.5, b = 1.5 in 90000 events for 5 and 
10 cascade steps. The wider curves correspond to 5 stages of the cascade. 
'Solid' curves represent the histograms with the correction p q (k) taken into 
account. 



